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Preface 


This study represents my effort to provide a method 
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for combining the responses of two modes of aircraft 


motion. The aircraft of interest was statically unstable 


aie 


which proved to be an interesting feature. The mechani- 
zation employed was demonstrated to have favorable re- 


sponse characteristics. 
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Abstract 


A root locus analysis and computer simulation are 
used to determine the feasibility of one proposed method 
of mechanizing the blending of the normal acceleration 
mode with the basic airerart response for the CCV YF-16. 
The root locus analysis predicts the stability and speed 
of response of the mechanized aircraft. The computer 
simulation confirms theseresults. Comparison is made of 
the responses of the basic, present CCV, and proposed 


mechanized YF-16. 


MECHANIZATION OF BLENDED Ay MODE 


FOR CCV YF-16 


I. Introduction 


Since the beginning of powered flight, man has sought 
to increase the effectiveness of his methods of control 
over his craft. Due to increased sophistication and demands 
for increased performance in fighter aircraft, improving 
aircraft control has become an important area of consider- 
ation. This has led, in recent years, to a greater emphasis 
on the use of more active controls and their implementation 
on Control Configured Vehicles (CCV). A CCV is one in which 
advanced control technology as well as aerodynamics, struc- 
tures, and propulsion are employed in the initial definition 
process (Ref 5:1). The Air Force Flight. Dynamics Laboratory 
and General Dynamics Corporation have been studying this 
concept on the YF-16 prototype airplane since December 1973. 
The CCV functions of direct lift and sideforce (constant 
angle of attack and sideslip, varying flight path), fuselage 
pointing (constant flight path angle, varying angle of attack 
and sideslip), and vertical and horizontal translation (vary- 
ing flight path angle at constant attitude) can be achieved 


with the use of an auxiliary flight control computer in an 


oe ee 


"add-on" configuration with the basic YF-16 flight control 
computer. The "add-on" configuration was necessitated by 
the program's requirement that the capability exist to re- 
vert back to the basic YF-16 flight control system at any 
desired point with no adverse flight characteristics. This 
also eliminated money and time costs in developing and cer- 
tifying a new control system. 

However, the “add-on" configuration caused a slight 
increase in pilot workload in certain flight regimes, e.g., 
air-to-air tracking, air-to-ground tracking, etc. Pilots 
are required to turn knobs and flip switches to transition 
from one CCV mode to another. The longitudinal CCV modes 
are available only through a two-axis force button mounted 
on the side stick controller. This non-optimum CCV button 
controller reduced the frequency of longitudinal CCV com- 
mands during flight testing. Years of conditioning plus 
neuromuscular considerations combined to make button inputs 
seem unnatural to the test pilots. Flight tests indicated 
that there was also a tendency for cross-talk between button 
and stick inputs. Alternate blending of CCV longitudinal 


inputs would alleviate this problem (Ref 9:19). 


Purpose 


The purpose of this study was to analyze one method of 


eee 


mechanization to provide blending of the CCV Ay (Direct Lift) 
mode with the control system of the basic YF-16 through side- 
stick force inputs. The mechanization was to provide a 


"floating" Ay mode. 


Scope 

The analysis in this study involved the blending in the 
longitudinal axis of the direct lift, Ay mode, to determine 
its feasibility and sy eames its stability. Two methods 
of analysis were employed for the purpose of this study. 
The first was a root locus analysis using the short period 
approximation of the YF-16 airframe combined with the flight 
control system of the basic YF-16 and the auxiliary flight 
control system of the CCV aircraft. The second method 
utilized a modified computer program simulation of the air- 
craft and its flight control system. The aircraft considered 
in both methods was the YF-16 prototype without the canards 
which were added to the actual CCV aircraft. This was done 
to simplify the analysis since the canards did not substan- 
tially affect the aircraft characteristics in the longitu- 
dinal mode. 

Evaluation of the results included both the root locus 
predictions and the time histories of the aircraft parameters 


for the basic YF-16 and the mechanized CCV YF-16. The 


: 
: 


evaluation covered two flight conditions; Mach = .6, alti- 
tude = 10,000 ft., and Mach = .8, altitude = 30,000 ft. All 
flight conditions were with the aircraft center of gravity 
at 35% MAC. Further description of the selected flight con- 
ditions is outlined in Table I. 

The physical description of the aircraft and diagrams 
of the CCV Flight Control System, including the Pitch Axis 


Functional Block Diagram are given in Appendix A. 


Table I 


Flight Condition Characteristics 


: 
; 
em 


Fat. Cond. 1 Fit. Cond. 2 

Altitude (ft.) 10,000 30,000 
4 Airspeed, Vp (Mach) -6 8 
(£t/sec) 646.8 795.88 
| Dynamic Pressure (1bs/£t7) 361.31 281.9 
| Air Density (slugs/f£t?) .001756 .00089 
Trim Load Factor (g's) 1.0 1.6 
? 

Trim Flight Path Angle (deg) 0 0 


Trim Angle of Attack (deg) 


II. Method of Analysis 


Assumptions 


The assumptions made in this study are as follows: 
1. The earth is flat and non-rotating. 
2. The aircraft is rigid and of constant mass. 


3. The atmosphere is at rest with respect to 
the earth. 


4. Gravity is a constant acceleration. 


5. The X-Z plane in body axes is a plane of 
symmetry. 


6. Any control deflections at equilibrium 
remain constant throughout the motions. 


7. For the root locus analysis, the short 


period approximation closely models 
the aircraft. 


Description of Root Locus Analysis 

The linearized representation of the equations of motion 
was used to obtain an indication of the instability of the 
basic airframe of the YF-16. From the linearized equations 
of motion, the short period approximation was formed in body 
axes in the manner described by Roskam (Ref 7:Chap 6). 

The linearized longitudinal equations of motion with 
dimensional derivatives employed were the X and Z force equa- 


tions and the pitching moment equation. These equations are 


as follows: 


ee 


1 


RE mn peer ~~ ne ocnwter A AEE AT 


ex, 


X equation - 
e ae . a . $ 
u Vpasin a) DEER Ht Sy Re Ree tae a 
+X, 5¢ (1) 
£ 
Z equation - 
VR" - VRdcos my = g@sin e, + 2 + Za + Zia + cs + 
+Z. 6 
£ 
: (2) 
M equation - 
= bi + a +M,a + Mp + Mo + Mgq + my 4 + me 


Equations 2 and 3 yield a short period approximation of the 


form: 
VRS ~ ah, -(V,cos a, + Zq)S -g oane a(S) A 
-(M;S + M,) s* - Mqs 6(S) 
Z 2 
. oe 
8 (s) + 5, (s) 
“ . ms 
e £ (4) 


The development of these equations through the short period 


approximation can be found in Appendix B. It is noted here 


pact a i st Ni tthe 


that the necessity for the inclusion of the flap parameters 


arises from the fact that flap deflection is the principal 


method by which direct lift is achieved for the Ay mode. 


It then becomes necessary to form the short period approxi- 
mation transfer functions with respect to flaps as part of 
the basic airframe dynamics. Equation 4 thus becomes the 
basic model of the aircraft dynamics for the root locus 
analysis. 

Since the YF-16, from which the CCV evolved, is a 
statically unstable aircraft, it, by necessity, is a highly 
augmented aircraft. To form the complete aircraft model 
requires the inclusion of the flight control system for the 
basic aircraft as well as the auxiliary flight control system 
for the CCV aircraft. Figure 1 shows the simplified longi- 
tudinal control system block diagram used in this analysis. 
The appropriate transfer functions are listed in Appendix C. 
The forward loop transfer functions G(s), G(s), and G, (s) 
are: 


G, (s) 


Go(s) = 


G3 (s) 


we 3sks [or}UOD AUSTIA 9T-dA OTSed peTJTIGuTS ~°T eanbTtJ 


G4, H Ho and H. are flight control augmentation transfer 


i" 3 
functions. These functions are formed either directly or 1 
indirectly from the short period approximation. Develop- 
ment of the a,(s)/ § ,(s) transfer function is shown in 
Appendix D. The development for a,(s)76 (8) is similar. 
The coefficients associated with all transfer functions in 
this study are combinations of scheduled gains already em- 
ployed in the flight control systems, the coefficients com- 


puted in the short period approximation, and the coefficients 


associated with the actuator dynamics. 


Reduction of block diagram gives a complete model of 


RE TERURT TEE 


basic YF-16 to be used in this study. Once reduced, the 
overall transfer function is a,(s)/g(s), normal acceleration 
to commanded g's. This transfer function should be stable 


and have acceptable characteristics at both flight conditions 


2 EE EA SEE RSET: 


considered. 

| Figure 2 shows the simplified block diagram for the 
CCV YF-16 in the Ay mode. At first, it appears difficult 

to analyze this block diagram. However, using block diagram 


algebra and Mason's Rule (Ref 2:151), the task is simplified 


Saidibe Seihinhicameiataihh dais seeecnie ation 


4 as is shown in Figure 3. The forward loop transfer functions 


G,(s), G,(s), and G(s) remain the same as for the basic 


NPN tee oT 


aircraft. However as the feed forward paths are closed, 


Perms nes, 


Re ahcilaneteto 


wo3sKs [TOr}UOD AYSTTA 9T-dA ADD pEeTFtTduts 
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weqzsXs [Or}UOD 3YUSTTA ADD UMeATpeu “°E€ eanbTy 


. 


a 


caution must be exercised to ensure that the denominator 
polynomials of the transfer functions with respect to 

5. (s) become those of the now augmented aircraft. As 
each feed back loop is closed behind the summing junction 
where the feed forward loop enters the main forward path, 
the aircraft essentially changes roots and becomes more 
augmented. The problem essentially reduces to a modifi- 


cation of the forward loop transfer function G as seen 


4’ 
in Figure 3. Using Mason's Rule and Figure 3, the a,(s)/ 


—§2(s) closed loop transfer function for the CCV aircraft 


becomes: 
an - G)G,G, (G+ Go + G+ G5)+ GG, (1+G,H,) 
ay = 
r) = + 
e 1 a G,S,H, G,6563H, 


+ Gg (1 + GH, + G GH.) 
(8) 


Equation 8 should indicate a stable CCV aircraft in the Ay 
mode. The transfer function and gains for the CCV aircraft 
are listed in Appendix E. 

The mechanization employed in this study is shown in 
Figure 4. The mechanization can be redrawn as shown in 
Pigure 5 to show how it functions. Between the stick force 


input and the Pitch Stick Gradient (shown in Figure 33 in 
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Appendix A) is placec a low pass filter which allows the 
low frequency commands to enter the basic aircraft flight | | 
| control system. Between the stick force input and the CCV 
Ay mode, the mechanization acts as a washout allowing the 
high frequency transients to enter the CCV path. As time 
increases, the signal is washed out from the CCV path 
leaving only the low frequency commands flowing into the 


basic aircraft. The gain, K associated with the mech- 


10’ 
anization will be determined from the flap deflection time 
history plots obtained from the computer simulation. The 
value of Kio will be determined by the response of the air- 
craft and the length of time that it is desirable to have 
the flaps deflected. The values tested are Kio = -l, 1.0, 
10 and 100. 

The limiter in the CCV path is required to ensure that 
flap commands do not exceed the travel allowance of the 


symmetrically deflected flaps. 


The an transfer functions, a (s)/g_(s) and a (s)7 8 Ws). 


can be used to obtain an indication of the aircrafts stability 
and response when mechanized. The an transfer functions can 

: | be combined to obtain the overall closed loop transfer func- 

: tion relating normal acceleration to stick force input as 


follows: 


SE ST 


a a a eee eee ang 


=) 
a 
Nt 


where 5 / 5 op is a gain that converts stick force input (in 


pounds) to a commanded flap deflection (in radians). 


Computer Program Analysis 


The computer used in this simulation was the Wright- 
Patterson AFB CDC 6600 computer. The program used for the 
simulation is basically that in Reference 1. Since that 
report was published, the program has been used in de- 
parture studies of the YF-16 by the Air Force Flight Dy- 
namics Laboratory. A listing of the program used in this 
study is provided in Appendix F. 

The main program provides for setting the initial con- 
ditions, calling subroutines to trim the aircraft (TRNTRIM 
and ANGLE) and to integrate the equations of motion with 
respect to time (subroutine RKGXYZ) and plotting the time 
history plots of important aircraft parameters. 

The subroutine RKGXYZ is a fourth order Runge-Kutta 
numerical integration scheme that integrates the equations 
of motion. The procedure computes 200 integration incre- 
ments per one second of time history displayed. 


The subroutine GYRATES is called by RKGXYZ and pro- 


a K a 6 a 
2) E-File 
SF S+Kio Ie sr} [S* Kio b (9) 


ct a Tt 


us 


vides the modelled flight control system, provides values 


of the aerodynamic coefficients using a table look-up 


scheme, and defines the equations of motion including the 


Euler relations. 


The equations of motion defined in the program are 


the full non-linear equations employing table look-up aero- 


dynamics. The equations are as follows: 


Equations of Motion 


2 
“ pvp S : 
U = -g sin@ + VR - WQ + oe c. * 208 ¢ 
pv, 
V=g cos @sin@® + WP - UR acy 
‘1 pv.'s T sin€ 
=: ea ioe ass 1 
W = g cos@cos @ + UQ VP + ans Co a 
1 £ =f ' 1-1 
p = ——t aaa | OF 4 == pol(i- a 
1 - —*42 q. I. I 
IyIz 
I 
xZ Sb 2 Sb - 
ae aca * (@ * V_-C 
5 2h? es i ie ae 
Zz x 
pe ie oy = ae R= p* ae Vp C 
v=) x ac. CUP Rk Se 
y y 
~~ * I - 2 
2 XZ Sb 
a Tr pea e B- OR|+ 57 pv,"C, 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


h = U sin@é - Vcos@ sing -W cos @ cos ¢ (16) 


Euler Relations 


© = Q cos @ - R sin ® (17) 

ie Eee ee Q sing + R coso| (18) 
cos © 

e e 

@’=P + © sin @ (19) 


Cc Cor Cie Cm, and C, are functions of angle of attack. 


x? Cy. z 

As mentioned earlier, the principal method of obtaining 
direct lift is through deflection of the symmetrical flaps. 
The computer program had to be modified to include the aero- 
dynamics of the symmetrically deflected flaps. Data was 
obtained in graphical form for the flaps at M= .2. This 
Mach number was chosen since the aerodynamics package with 
the program contained data at M= .2. That is, any flight 
condition for which time histories are computed uses aero- 
_@ynamics at M= .2. Therefore, any errors due to change in 
Mach number would at least be consistent. Given time his- 
tories that were close to actual flight traces for the 
basic YF-16, it was found that the errors at increased Mach 
numbers were acceptable. 


The symmetrical flap aerodynamic data was obtained from 


Reference 3, Page 18, in graphical form. The change in lift 
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coefficient due to flap deflection was found to be essen- | 


tially linear in angle of attack between a = -4° and 
a= 12°. The function form of this change is as follows: 
AC. = .016 6 (20) 


L £ 


The change in drag coefficient due to flap deflection was 
found to be composed to essentially two linear segments. 


The functional forms are as follows: 


AC, 
= ,0002 + .0001958 a a 20 (21) 

A 5, 

Acy 
= .0002 - .000325 a a <0 (22) 

7A) b¢ 


The change in pitching moment due to flap deflection was 
found to be of the form: 


Ac 
pera 
Aéb¢ 


= (.0027 - .000659 5 ¢) + [<--00748 + .0000042 5) a] 


[(--000052 + .0000277 5.) a*| + [‘-000012 - .0000015 5.) oa 


- .000000403 a (23) 


These functions were found using a least squares curve fit 
(Ref 4:Chap 7). 


During the simulation, time history data for the basic 
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YF-16, CCV YF-16 and the mechanized YF-16 at the two chosen 
flight conditions were plotted together to aid in evaluating 
the various responses. For the mechanized aircraft, the 
various values of Kj) (Fig 4) were inserted to compare the 
lengths of time that the flaps were deflected. The most 
reasonable value was then used in the simulation for the 


mechanized aircraft. 
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III. Discussion c 


Root Locus Analysis 

The non-dimensional stability derivatives were deter- 
mined for the two flight conditions and are listed in Table 
II. Using the definitions of the longitudinal dimensional 
stability derivatives found in Roskam (Ref 6:6-17), Table 
III was prepared. Table III and equation 4 lead to the 
basic airframe short period approximation transfer functions 
for each flight condition. These transfer functions are 
listed in Appendix E. For the 4, transfer functions, the 
Root Locus plots are shown in Figures 6 through 11 for 
unity feedback around each transfer function. The instability 
of the basic airframe is clearly observed in all plots. 

The above transfer functions were combined to obtain 
the forward loop transfer functions G,(s), Go(s), and G3(s). 
These coupled with the control system outlined in Figure l 
formed the model of the basic aircraft. The root locus plots, 
as each loop of the control system is closed, for both flight 
conditions are shown in Figures 12 through 17. It can be 
observed that simply closing the pitch rate feedback loop 
results in a stable augmented aircraft. The closure of the 
remaining two loops simply improves the performance without 


impairing the systems stability. Closing the last loop and 


Table II 


Stability Derivatives 
(Nondimensional, stability axis system) 


Flt. Cond. 1 Flt. Cond. 2 


Cc -005 -020 


c «£95 - 280 


Table III 


Stability Derivatives 
} (Dimensional, Body Axis System) 


eee UI EEE nnn SEI RINEER IES RERESREnenemnennnmemmnmtnel 


Flt. Cond. 1 Flt. Cond. 2 
eee Sel, fap hee as Sno a a ks ce ae ee. oh Pe Rh 
Z, (ft sec 7) -721.8 -600.5 
Zo (£t sec ~+) -1.895 =3.45 
Zz (ft sec 72) ~103.7 =68.-22 
b, 
| an 
Zs (ft sec ) -125.3 -~106.03 
£ 
M, (sec -2) 4.02 1.70 
| M, (sec ~+) =~.327 “= 081 
j 
Mq (sec ~+) ~.743 -.523 
; M, (sec -2) ~j4.23 =11.16 
e 
M, (sec ~2) #435 =1.86 
£ 


cascading with G, yields the Normal Acceleration to G com- 
a 
manded closed loop transfer function, Bae)" listed below 


for both flight conditions. 
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Flight Condition Il: 


a, (s) (6.33) (s +5) (s +1) (s +10) (s ¥15) 
ee (s + 8.3) (s +1.368) (s + .6434) (s + 2.765 =3.389}3) 


i 
; 
= a 


(ea 2.021. 21.684) 
(s+10.7) (s 12.71) (S + 16.11 =15.89) (24) 


Flight Condition 2: 


a's) (7.026) (s +5) (s +1) (s +10) (s #15) 
Posius (s + 8.3) (s + 10.55) (s * 12.83) (s + .7696 ~.3564}) 


SS NRA 


i (s +1.21 = 9.4643) 
(s + 2.643 43.323) (s + 16.13 +15.133) (25) 


The various gains (as defined in Appendix C) are listed in 
Table IV. 

From Figure 3, it can be seen that the CCV aircraft 
would not present any new information from a root locus 
point of view since the basic aircraft portion remains in- 
tact. The forward transfer function Gy is essentially what 
is being modified in this configuration. Therefore, the 
Normal Acceleration to Flap commanded transfer function, 
is), was found directly using Mason's Rule. The appropriate 
values of the gains are also listed in Table IV. The transfer 


functions are presented below. 


op 


RF A NE, 


A ree UE OT 


ee 


Table IV 


Transfer Function Gains 


Flt. Cond. 1 Flt. Cond. 2 
GK, 284.56 223.2 
GK, 0112 .0076 
GK, »265 479 
K, (basic) ~7.46 -8.54 
K, (CCV) ~.525 -.473 
2 
K. .758 .665 
Ke ~.426 -.701 
K, -.362 -.658 
Ke -15.52 -10.64 
Kg 4.72 3.58 
Ky 1.079 1.235 
Ky 5 5 
K 40.46 46.29 


Flight Condition 1: 


a, (s) (.4706) (s + 3.18) (s +1.231) (s + 9.978) (s+ 15) 
be (s + 3.229) (s+ 1.367) (s + 3.987) (s + .642) (s+ 12.72) 


s +15) (s+7.557) (s + 24.17) (s + .8539 +.86283 
(s + 10.7) (s+ 8:3) (s + 3.144 + .024833) (s + 2.767 + 3.3833) 


(s+ 1.004 + 11.683) 
(26) 


(a +16.11 £ 15.974) 


Flight Condition 2: 


ay (s) (.9499) (s + .7417) (s + 10.02) (s + 7.849) (s + 15) (s + 15) 


S, (s + 3.984) (s + .769) (s+ 10.54) (s +15.09) (s +12.84) 


(s + 18.88) (s + .699 + .60173) (s+ 1.206 + 9.4563) 


(s+ 8.3) (Ss +2767 * 3553) (8 + 2.653 = 3.3159) (eo + 16.12 


+16.16}) (27) 


For the mechanized aircraft, the two transfer functions 
(Eqs. 24 and 26 for flight condition 1 and Eqs. 25 and 27 
for flight condition 2) were blended according to Figure 5. 
The closed loop transfer function was formed for values of 
Kio = 0.1, 1.0, 10.0 and 100.0. This transfer function re- 
presents the response of normal acceleration to stick force 


inputs. Figures 18 through 25 show the time response for a 
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unit step input for the various values of K 


10 


flight conditions. The settling time of the various re- 
sponses are listed in Table V. It can be seen that the 
higher the value of Kyo: the faster the settling time. 
However, from Eq. 9, it is determined that the higher values 
of Kyo are in effect causing the basic aircraft response 

to dominate and penalizes the Ay mode response by washing 
out the flap command at an increasing rate for increasing 
values of Kio: Therefore, it is desirable to use the lower 


values of Kio: This will slow down the response of the 
basic aircraft, increase the time the flaps are deflected, 
and allow the flaps to reach a high maximum value which is 


closer to the actual CCV flap deflection. This will be 


shown more clearly in the results of the computer simulation. 


Computer Simulation 

As mentioned earlier, the characteristics of the CCV 
Ay mode were faster response in normal acceleration and an 
almost constant angle of attack while rotating the velocity 
vector. Therefore, in simulating this mode, it was apparent 
that the parameters of importance were angle of attack, 
normal acceleration, and altitude as a function of down 
range distance. These parameters can be compared to eval- 


uate the response of the various modes. 


Table V 


a 
Settling Time (sec) for _n_(8) gor values of K 


sf +0 
Kio Flight Flight 
Condition Condition 
1 2 
0.1 40.66 40.47 
1.0 7.24 5.80 
10.0 5.76 4.04 
100.0 5.67 3.94 


The inputs for each run of the program were trim angle 
of attack, initial altitude, trim load factor, and a history 
of stick force inputs at discrete time intervals. The pro- 
gram would then give the trim conditions and initialize all 
the state variables before it worked its way through the 
flight trajectory. The scheme for each run was to input the 
trim conditions specified in Table I and after one second 
introduce a positive six (6) pound step input. Data for 
normal acceleration, angle of attack, and altitude was col- 
lected for plotting of the basic, CCV, and mechanized air- 


craft responses at each flight condition. 


Sea $$ $e 


Once data was obtained for the basic and CCV aircrafts, 
the two modes were mechanized according to Figure 4. The 
value of Kio was allowed to vary in successive runs and 
the time history of the flap deflection was recorded. For 
a particular flight condition, the time histories of flap 
deflection far the various values of Kio were plotted to- 
gether with the flap deflections during just the CCV mode. 
These plots are shown in Figure 26. It is clearly obvious 
that the lower value of Kio (0.1) yields the greatest de- 
flection in flaps and, consequently, the more pronounced 
response of the Ay mode. The time over which some portion 
of the flaps remains deflected is increased for Kyo = 0-1. 
This is in accordance with the prediction of the root locus 
analysis. 

The important parameters listed at the beginning of 
this section were then plotted together for each flight con- 
dition. The value of Kj) chosen for the mechanized aircraft 
was 0.1. The plots are shown in Figures 27 through 32. It 
is observed that the mechanization of Figure 4 allows a much 
quicker initial response and then evolves into the response 


of the basic aircraft. 
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IV. Conclusions and Recommendations 


Conclusions 

The analyses of this study have led to the conclusion 
that mechanization of the blending of the CCV Ay mode and 
the basic YF-16 flight control systems is both feasible 
and worthwhile. The mechanized aircraft has effectively 
combined both modes to give increased performance over the 
basic aircraft. The mechanized function is available to 
the pilot through the sidestick controller and does not re- 
quire inputs from a two axis force button. The cross talk 
between button and stick are eliminated. If the command is 
held as a step input, the CCV Ay mode will washout leaving 
only the basic aircraft mode. If faster inputs are required, 
the CCV Ay mode will predominate and little or no response 
will be seen from the basic aircraft. This is clearly an 


advantageous system for aerial maneuvers. 


Recommendations 

This study has only touched upon one possible method 
of blending the Ay mode with the basic aircraft. Other 
untested methods have been proposed by some of the aircraft 
industries. These methods should be tested and compared 


to the method of this study in an attempt to find the most 


practical solution. 

Five more modes are available on the present CCV air- 
craft, of which all the longitudinal modes are available 
only through the two axis force button. These, too, could 
be blended in some fashion. 

There are almost limitless possibilities of improving 
the performance of modern fighter aircraft. To maintain 
the advantage of air superiority, these possibilities must 


be explored. 
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Physical Description of the Aircraft 


602.5 slugs 


= 260 ft? 


29.0 £t 


10.94 ft 


8400.9 slug - ft? 


47000.0 slug - £t2 


55000.0 slug - £t* 


550.9 slug - £t? 


Control Surface Deflection Limits 


Biavetor . <-a.+ s% «= « 22S" 


Plapevons «444 « % 220° 


Rudder .96 2 @ = 6 < « & +30° 


: 
; 
: 


WERE TO 


COV PITCH AXIS AND 
CATERAL -QIRECTIONAL AXIS. 
FUNCTIONAL BLOCK DLAGRAMS 


ewe ene en ene eee 


ESTO Fonct EHOW 16 16) ' 
‘ 
Secor (DEG! FROM YF-16 FCC 
' 
' 
1 
'@ BIAS FROM @) 6 @2 
MODES TO YF-10 FCC ! 
' 
1 
‘ 
i] 


BA, (G's) FROM YF-16 FCC 


‘STICK FORCE FROM cc 


[Yaw TRIM (EG &y ) FROM vE-16 FCC) 


RUDDER PEDAL FORCE (LB) 


FROM YF-16 FCC 


Wb (oecistc? From vF-16 FCC 


SIGMALS FOR CCV 
SMASTER SWITCH DISCOMMECT 


‘hy OAS TO YF-16 FCC 


ce (Gs) FROM YF-16 FOC 


ere men oe ewe ow mew een wn ewe we oo ee ne er enn ore we we = ee 8 et er er enn ee wen me eH oe fe ee ee ewe eee ee eres errr eres — 


| 
j 


r 


THIS PAGE IS BEST QUALITY PRACTICABLE: 


wr wm a rn oo rn en ne one eee: 
Pa eereerese: 

SIGNALS FROM STICK COMMAND PATH OF YF-10 FCC PITCH STICK C COMMAND 

TO CCV Ay MODE S OVERLOAD PROTECTION CIRCUITRY | FROM YE-16 FCC PRIOR TO 


wc oe REDUCTION BY STRUCTUMAL 
es STICK G BIAS FROM Ay, & @2 MODES 10 on ieee 


COMMANC 16 4) PITCH STICK COMMAND PATW OF YF -10 FCC 


PITCH TRIM 
COWMAND 


Att O.00y EC) 


PITCH STICK FORCE 


TRUE ANGLE-OF-ATTACK 
@, 16 


@ mom 


PITCH Rare aa “ 60) 


mite CCW MASTER SwiTCH ENGAGED 


A) MLICTION OF any *UBai COmmnnonD 
CC¥ MODE HITCHES OUT MOtenn, POR 
COmman Pat 10 eUbees 


i wre 13 1G MOC) O toe 7) 
: — zi 
ROLL STICK FORCE 7 - + efie) al 
et ES oe 
t u-- ---~]- stenc. aang c---4 | 
tO aa =a Ditone st es] +48 
‘ r) (EGA) steucToNaAn tu ne % 


' YAW RATE 
: WY vote we, 
wom tou. a8) G) 


LATERAL ACCELERATION 
ow 


See eee ee ne ee ey 


Lwaencwmenon 


fiIS PAGE IS BEST QUALITY PRACTICABLY 
FROM COPY FURNISHED TODDG =__ 


en nn ne oe we eee ee wa ee oo a i ao ee ee ee ee a aw ay 


cey 


patton ta 
YF-16 FORWARD LOOP INTEGRATOR 1 & F2 SCHEDULED CAINS USED IN YF-I6 


| “ 
"by /ac FCC FOR STRUCTURAL OVERLOAD PROTECTION 
OS*ag/ac I 
SIGMAA FROM YF-10 FCC USED TO REDUCE | | oe | FROM PITCH AUTOPILOT, @2 MODE, 


Th PITCH STICK CG COMMAND AND TO OR A, MODE PURIFICATION 
PROVIDE STRUCTURAL OVERLOAD PROTECTION 4, COMMAND TO YI 


ics 
BY STRUCTURAL 
OVERLOAD PROTECTION SIGMAL 


FE SCHEDULED GAIN USED 
1M Va AXIS OF YE-ie FCC 


‘ 
‘ 
‘ 
' 
' 
1 

J 


on we ee ww ee ow ww | pe ww wm oo ee re ew wn oo wo | oe we we ow = oe | 0 TS ORES OS Hmmm 


PITCH AXIS 


STATIC STABKLITY 
omy 


NORMAL RUDDER PEDAL COMMAND SWITCHING LOGIC 
Py STATIC manssuat FOR Gran merssuan ALTITUDE 
© STATIC PRESSURE AT SEA LEVEL (29.9713" ny STO 
= 10tay PanSiuae 
= manact PRESSE OF, ~ 
A = ROH NueeeE, e_ 


Joie CC¥ MASTER SmITCH ENGAGED FOLLOWING OIENGAGEMENT OF CC v MASTER BeITCR 


MAMUALLY OB BY AUTOMATIC OSCONNECT 


Ay MUECTION OF AMY PiRay COmmaNOEO 
CC¥ MOE SHHTCHES OM reDMnna AEOAL | a: MOMMAL FIDAL COMMAND RimAINS SWITCHED OF 
Command Mate 'O meKe UNTR PIDAL FORCE 1) MOUTED SUCH Tmat Tot 

QUOT COmmAND 15 ESS Tran 2 DEGREES. rats 


«hs CC ROtARE PRGAL MLECTOR KHOR 15 
TuanaO TO Trt DBF POSITION. OF 


2. CC¥ OFSLONMECT BUTTON (wore MEL) One 
SIC GOH 15 DEPREIHO AND MELO 


ts? + TS + 52) ‘ ‘ be, 
a) ' RH HORIZ TAIL 
[= -t3* 


A. 


D 
it d /| UM WORIZ TAIL 
ut ‘ a 


30 OIG WC 


{Mart naet 170 F674 | 


MuUCTOR 


Frere eee se 


RH CAMARO 
SERVO COMMAND 


REFER TO 


CCV LATERAL-DIRECTIONAL AXIS: 
FUNCTIONAL BLOCK DIAGRAM 


* OQ) se ug? 
S ae ; 
BSED SRE Cipey ON Rare ety are eee Belted wa hited eat 


-|----------}--- = }-- - ed 


‘Ut CAMARO 
SERVO COMMAND 


3, COMMAND TO YF ie fcc 


1 
' 
7 {UR* f.28 COMMMAND TO _YF-16 FCC 
' 
' 
' 


ee ee ne 


22 gh 


rs THIS PAGE IS BEST QUALITY PRACTICABLE. 
FROM COPY FURNISHED TODDC ___— 


eee 
FL GFZ SCHEOULEC GAINS USED im YF-10 
| FCC FOR STRUCTURAL SvEMLOAD PROTECTION 


J 


FU SCHEDULED GAIN USEO 
1m YAM AXIS OF YF-10 FCC 


‘ 
! 
‘ 
' 
' 
t 
J 


> 2 © © © \0© 18 « «se 19 1 oe 
| boienqe 1 ' ‘ 1 Q t poste. i t + | 
eee SS ee ee ee 


H Uc ~ese os 
e s Eee Se ee eee e 4e cits it : : ' \ Bn iets { 
pat as (a a 
| - } i ' t i ; 
TANNOY eAI}--——----.-- ; = Reems nd 
By = Sta meassuan (08 Given pmasiuan ALTITUOE r H i ‘ i ' ‘ 
+ SANG pra ssumd At MA EVEL (20.0919" iy S10) ‘ = 0 at 1206 pol ote AG tardies DOWN H i H ; 
+ x lb seas ' Qe FST vise ow canter ‘ jp | om nn. a 
& 
” 


us 1s we 
-s PITCH STICK FORCE 


rote 
\/ \~u tose fe 
bn, #08 TRO HOMO; Tame) + 
ee womiz tant FIG. A PITCH COMMAND GRADIENT ipa a 
we “4 (CLUES BMAROUT FORCE) 


ee xusnocnismente 
PICLTER) Be” FLIGHT CONTROL SYSTEM FUNCTIONAL BLOCK DIAGRAM 
iia a en a Figure 33. (Ref 6) 
62 = 


THIS PAGE IS BEST QUALITY PRACTICABLE 
FROM COPY FURNISHED TODDC ___— 


PITCH TRIN COMMAND 
To rf-te FEC 
Lo 1ROM STICK COMMAND 
latin Pua tatu OF YF-ie SCC TO CCY 
ay & @2 noms MOOL_C OVERLOAD FROTE! 
TO PITCH STICK ComMAaD Checcttay @sias from aise. 
PATH OF WE ie FEC —eC moees TO VF-ie FCC 


pitce Tou rivce sticn 
comme (c's) CBOnme (C's) 


PIICH STICk FORCE 
PHOM Vr-Le FEC 


®ooy (te) 
FROM Ray COMPUTER 


sent LE NE 340 
roe 


| 


Ses 


nate < es 


cow masta 
wee 


Piree sr 1ce ¢ comae 
Fea 1418 FCC Perce TO 


Cae 
BLN. PREETI Sora 


ui 


't(iiiierp_m.. piciateh 


WHS Paez Is BEST QUALITY PRACTICABLE 
FROM COPY FURNISHED TO DD. 


C7) CC¥Y = AUTOPILOT 


DIRECTIONAL 


ATMOSPHERIC DATA 
FROM YF~1b FCC 


SIOMALS FAN STICK Conmano (owe urrecasroe 7 . +| 
PaTW OF YF ie ca 
MOOR C OVERLOAD PROTECTION fem ertce aurorivot. @2 moos, 
cla iTey @uias room ais a. 
MORES To VF-1e FCC 08 yy PORE PURIFICATION 


CoY PiTt 


APFOL / GENERAL OCYNAMICS ? 
Figul 


i IS BEST QUALITY PRACTICABLE 
% FURNISHED TO DD 


CCV = AUTOPLOT | LONG | OIRECTIONAL ® 


™® ae PUR 
tae 
Or '@ 


F 
'® a 
e RNOR]| MAN ENH 


u 
vi = Gd tie 
@2.A,, NODE @i moDE + 2 
ATMOSPHERIC DATA ag,  , Saenteeee : i ‘ ay, i ae, 
FROM YF-1b FCC ys dee pie 


— te +t iy 
Se ee 


Se ae 


- 


121075: CCB/RRG 
012776: wCw 


031076:CCB/RRG/WCW 


UOT Leek CCY PITCH AXIS FUNCTIONAL BLOCK DIAGRAM 


stCa 
APFOL / GENERAL DYNAMICS 


Figure 34. (Ref 6) 


63 THIS PAGE IS BEST QUALITY PRACTICABIEE 
FROM COPY FURNISHED TODDC ___- 


ae 


Appendix B 


Starting with the general non-linear equations for the 
forces and moments in component form 


i component of force: 


nf +QW - ve | = mg_* 


component of force: 


n|v+ up - ew] = mgy +*Fa 


component of force: 


n|@ + Pv = 
component of moment: 
XZ 


= tF° + RQ(I, - = 


j component of moment: 


5 2 2 
Qr, + (1, - 1,)PR “io (P -R ) 


component of moment: 
= PIL. + (I, - T,)PQ +I OR = N,+N,, 
The angular velocity vector is defined as: 


@ = Pi +Qj +Rk 


Using Euler angle Eq.(B7) becomes 


®= i(-Wsine + @ )+j( cos @singe . cos #) 


+k(W¥cos@ cos - 6 sin¢) (B8) 


Equating the components of Eqs.(B7) and (B8) yields 


P=6-wW sino : (B9) 
Q = Wcos@sinée + Ocos® (B10) 
a ; : 
; R= ¥ cos @cos®@ - O@sin& (B11) 


| The gravity vector is defined as 


9 = g,i+9 5 tazk (B12) 


Using Euler angles Eq. (B12) becomes 


g = i(-g sin®) + j(g sinécose) »K(g cos @cos @ ) (B13) 


Equating the components of Eqs(Bl2) and (B13) yields 


o, * 9 sind (B14) 
Jy = g sinécose@ (B15) 
: = g cos*#cos Oe (B16) 


Eqs (B14),(B15) and (B16) can now be substituted into 


Eqs(Bl), (B2), and (B3) 


| m(U - VR*QW) =-mg sin@+F, +F B17 

| ( QW) g eg T. ( ) 

| m(V + UR-WP) = mg sin@cos@ +F, +F (B18) 

B Ay Ty 

| e m(W - UQ + VP) = mg cos @cos6 + Fa. * F (B19) 
zZ Tz 


Using the perturbed quantities 


U=U,*u V=V,*F 


1 
P PL*P Q=Q,+4 


} + 6 
% 8 o,* 


assuming the peoeiates quantities are small, assuming small 
angles, eliminating the steady state terms, and imposing the 
restrictions 
pe Se hy n= eae 

the linearized equations are formed. Using the reldtions 
in body axes for the following 

= Vp cos a] 

aN. sin a, (B21) 
The equations of motion for the longitudinal axis become 

X force: 


m u+Vp aqsina, = -mg @cos 6, + 


Z force: 


m VR a - Vp dcos a, -mg @sin@, ar rs fr. (B23) 


M moment: 


Iv q= — mn 


Rewriting Eqs(B22), (B23), and (B24) with dimensional deriva- 
tives in the force terms yields 
X force: 


a= i - 6 + 
u Vad sin oF g eo e xu ax Ske 


u 
: X Se ox. Sf (B25) 
i e £ 
.| Z force: 
Va%- VRq cosa, = g946sin 8,* Zju+Zya * Zod 
+ “4 oi ae Re: (B26) 


M moment: 


Ge= Mu a ae +M o*My a@ +M.oa + M 
qa 


Taking the Laplace Transforms of Eqs(B25), (B26), and (B27), 


collecting terms, and writing in matrix form yields 


(ay - Fy? “Ks (Vsin a,s +g cos @,) u(s) 
-Zy VRS - 2, -(Vpcosa, * Z,)s-gsin9, a(s) 
(-M,-Mp,) — - (MgS *Mq ) s* - MS @(s) 
x 4 
‘ bp 
6 
Z, 6, (s) + |Z, (8) (B28) 
e £ 
M M 
be be 


67 


Taking the short period approximation (u(s) = 0) yields the 


result 


ee vps ~ De - (V,cos a, *ZQ)s -g sin6, « (s) E, 


2 
- (MS =M Se sad 6 (s) 


Appendix C 


For figures 1 and 3, the transfer functions have the form 
GK , 


1 
G(s) = 3350 Bish (cl) 


e 


G(s) = GK, yee (c2) 


& -(s) 


3 als (€3) 


G3 (s) GK 
K4 (s +5)? 


S (s +8.3) (ca) 


G, (s) 
Ks (s +5)? 


@5(8) = Tesi) (e <4) (e © 15) (cs) 


6 
G(s) (C6) 


G(s) (C7) 


G, (s) (C8) 


Kg 
(s + 4) (s +20) 
K, (s +5)? 


CES ICES) is 


Gy (s) (C9) 


H, (s) —_—— (C11) 


Hy(8) = 505415)? se 


The gains are defined as follows 


GK, = 20 
L “as 


GK, = avs (c14) 
e 
Cans; 
GK3 = Cy, (C15) 
e 
K, = (F1F2 - 1) (8.3) (1-25) (F3) for Basic YF-16(C16) 
G 
= (FIF2 - 1) (8.3) (1-25) (F3) for CCV YF-16 (C17) 
| Ke = 6 F3 ae x (C18) 
} eee be 
en Oe (C19) 
Cay, 
K, = 4 £ (20) 
q ae 
; e 
| Ky = 80 C (21) 
| Yb— 
| K, = 80 C (C22) 
: ake 
K, = (1.5) (F3) (C23) 
K, = (10) (F4) (C24) 
K, = (56.25) (F3) (c25) 


where Fl, F2, F3, and F4 are scheduled gains that depend 
upon the flight condition. 


A "Cc" indicates the coefficient associated with the 


' transfer function specified by the subscript. 


3. 


Appendix D 


2 SRR ae oe Ne RATTLE 


Referring to Figure 35, for $= 0, 


6= ary (D1) 


= 
a 
) 
8 
o 
R 
i} 


angle of attack 


and Y= flight path angle 
Therefore, 

| Y= 6-a (D2) 
Also, Vo=Vv a, (D3) 
and, vV=%¥ ea é, (D4) 
Since ey = cosy¥ i+sinyj (D5) 

implies é, = -ysiny i + Y¥cos7j = ¥(-sinyi+cos7j ) 
(D6) 
However, en = -sinyi+cos7j (D7) 
therefore, ey = ¥ en (D8) 


Substituting Eq(D8) in Eq(D4) yields 


Figure 35. Reference System for Determining Normal Acceleration 


<te 
" 


Ve, + V ye, (D9) 


Linearizing Eq(D9) in the normal direction yields 


Vv o= (U) +u) ( 7) +7) (D10) 
However, un = 0 and u is negligable 
leaving Vyr= u,7 = Ay (D11) 


normal acceleration about the center of gravity. 


From Eq(D2), ; ; 
Y= @=4 (D12) 


Substituting Eq(D12) in Eq(D1ll) yields 


A =U,(6@ - 4 D13 
we a) (D13) 


Normal acceleration at the pilot's station (location 
of the accelerometer) becomes 

A =A sia a (D14) 

Mec. “co, = 

where q = 6. 
Substituting Eq(D13) in Eq(D14) yields 

AND s. = 5, ( q-—-a@ ) ei 4 (D15) 
Taking the Laplace Transform of Eq(D15) yields 

a, (s) = U,( a(s) - s a(s)) +18 q(s) (D16) 
Substituting the transfer functionsq,, (s) and a>. (s) 
e 


e 
and regrouping terms leads to the result 


*n ig) = (1,8 *U,) as) - 0,8 a(s) (D17) 
3 ; 


e 


The accelerometer output equals the total acceleration minus 


gravity. However, the perturbation contributions due to 


gravity fall out at 0 = 0, ® = 0. 


Appendix E 


Short Period Approximation 
Transfer Functions for 
Basic Airframe 
Flt. Cond. 1 


— (s) = -.160 (s +89.48) 
5 (s- 1.0) (s +3.18) 


Ts) = = 14.228(s *1.17) 
(s- 1.0) (s + 3.18) 


Ang) = 720424(s*+ 2.022 5 +137.53) 
os (s-1.0)(s 2.09) 


a.) = =+194(s +7.68) 
(s) fe - 1.0) (4 © 3.16) 


& _ -=-1.287(s +1.78) 
c=) (s - 1.0) (s + 3.18) 


Bnis) = +059(s* + .935s - 13.58) 
8 (s - 1.0) (s +2.09) 


Short Period Approximation 
Transfer Functions for 
Basic Airframe 
Flt. Cond. 2 


~.085(s + 130.2) 
(s - .62) (s + 2.09) 


-11.16(s *.769) 
(s - .62) (s + 2.09) 


-.0407(s2 + 2.42 + 91.03) 
(s - .62) (s +2.09) 


-.133(s +14.4) 
(s - .62) (s + 2.09) 


-1.835(s 4.888) 
(s - .62) (s + 2.09) 


.0447(s* + .566s ~ 15.74) 
(s - .62) (s + 2.09) 
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Vita 


Captain Kenneth Race was been th Dover-Foxcroft, 
Maine on August 17, 1950. He graduated from Peru Central 
Jr-Sr High School, Peru, New York in 1968. He received 
a Bachelor of Science degree in Aeronautical Engineering 
from Rensselaer Polytechnic Institute, Troy, New York in 
1972. He received his commission through ROTC. His 
initial assignment was to the 32lst Strategic Missile 
Wing, Grand Forks AFB, North Dakota as a Deputy Missile 
Combat Crew Commander. Before his tour ended, he served 
as a Standboard Evaluator and a Combat Crew Commander. 

H He received his Regular appointment in 1975. In 1976, 
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